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Abstract 

In this Letter, we show that a three-dimensional Bose-Einstein solitary wave can become stable if the dispersion law is 
changed from quadratic to quartic. We suggest a way to realize the quartic dispersion, using shaken optical lattices. 
Estimates show that the resulting solitary waves can occupy as little as ~ 1/20-th of the Brillouin zone in each of the 
three directions and contain as many as N = 10^ atoms, thus representing a fully mobile macroscopic three-dimensional 
object. 
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1. Introduction 

Creating mobile self-supporting three-dimensional (3D) 
matter waves—3D analogs of the ID soliton^ first realized 
in Refs. [a [3] —is a long-standing goal of physics of ultra¬ 
cold gases laiiiiiiziii]- 3D Bose-Einstein solitary waves 
in continuum space are unstable because they violate [7] 
the Vakhitov-Kolokolov (VK) stability criterion 1^. One 
strategy for circumventing this problem is to create dis¬ 
crete solitary waves in optical lattices. These objects may 
be stable for some sets of parameters mala, but are lo¬ 
calized on a limited number of sites. Thus they occupy a 
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solitary wave is an isolated wave that maintains its shape due 

to a balancing of dispersion and nonlinear attraction; in strict usage 
[J, a soliton is a solitary wave with a special property: when it 
collides with another local disturbance—e.g. another soliton—then 
asymptotically far from the collision, it regains its initial shape and 
velocity. Regardless of how they may be realized, 3D solitary waves 
are not expected to be solitons, because they are not expected to 
instantiate any integrable systems. 
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substantial portion of the Brillouin zone and so have lim¬ 
ited mobility. Another stabilizing strategy is to oscillate 
in time the interatomic coupling constant around a nega¬ 
tive value, but switching to positive for periods of time [B] 
(the VK criterion does not apply to time-dependent non- 
linearities). However, here substantial atomic losses will 
limit the lifetime of the 3D object. In a third approach, 
3D solitary waves are stabilized by varying the nonlinear¬ 
ity strength in space [5]; but then they are not free objects 
moving in a translationally invariant medium. 

Here we propose a way to satisfy the VK criterion by 
changing the atomic dispersion law from quadratic to quar¬ 
tic. We suggest using shaking to couple the lowest and the 
first-excited energy bands in a 3D optical lattice so that the 
quadratic portions of the dispersions cancel out, and the 
quartic terms become dominant. This results in a highly 
mobile localized object. 
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2. Dispersion law and the stability of solitary mat¬ 
ter waves 


Let us discuss how the dispersion law enters the VK 


criterion. Our main conclusion, given in Eq. (14) below, 
may also be reached using the well-known heuristic argu¬ 
ment based on how the kinetic and interaction energies 
scale with the typical length (we will summarize this rea¬ 
soning following that equation); but because of the nov¬ 
elty of the situation, we prefer to use the more rigorous 
route via the VK criterion. Consider a time-dependent 
nonlinear Schrodinger equation for the wavefunction of a 
d-dimensional Bose-Einstein condensate with a generalized 
kinetic energy represented by a particular differential oper¬ 
ator of p-th (p being even) degree (below we will propose 
a scheme, using ultracold gases, for realizing a physical 
system described by Eq. ([^ with d = 3 and p = 4): 

d 




dt 


( 1 ) 


i=l 

where p > 0 is the dispersion strength, p < 0 is the cou¬ 
pling constant, and pi = —ih-^ is the operator of the ith 
component of momentum. Assume that Eq. ([^ admits 
normalizable stationary solutions 


4'(r, t) = exp[—ipt/h]tj;{r) : (2) 


their spatial part ip{r) will then be governed by a time- 
independent nonlinear Schrodinger equation, 

d 


p > 0, p < 0. 


( 3 ) 


Any solution of ([^ is a point of extremum, subject to the 
normalization constraint that 



( 4 ) 


be the number of particles, of the energy functional 


Here 


ap/2^ 

r = )«!■i'r Ht 


d „ 

T = r]hPY^ j(Yr 


and 


( 6 ) 

( 7 ) 


are the kinetic energy and the interatomic interaction en¬ 
ergy functionals, respectively. The corresponding varia¬ 
tional space is the space of continuous functions of the 
coordinates; the chemical potential p enters as the La¬ 
grange multiplier enforcing the constraint in Eq. We 
can conjecture a stability criterion for localized station¬ 
ary solutions of Eq. ([^ already at the level of the virial 
theorem: regarding the interaction functional in Eq. ([^ 
as a mean-field approximation of the pairwise interaction 
energy with the d-dimensional d-function as the potential, 
the virial theorem m would predict pT = — dV, and sub¬ 
sequently 


E = 



( 8 ) 


Now, one would expect the energy of a stationary state to 
be below the energy of a dilute, heavily delocalized cloud; 
the energy of the latter is close to zero, so the right-hand- 
side of Eq. (|^ must be negative. At this point we con¬ 
jecture that for solitary waves to exist, the dispersion law 
must be sufficiently sharp, namely p > d. Next, we will 
confirm this using the rigorous VK criterion. 

First, observe that if one multiplies the stationary non¬ 
linear wave equation in Eq. (|^ from the left by 'L*, inte¬ 
grates it over all space, and combines the result with the 
expression for the energy given in Eqs. ([^|^[^, one ob¬ 
tains [7] E = fiN — V. Combining this with Eq. ([^, we 
obtain a relationship between the chemical potential and 
the kinetic energy that, in particular, provides information 
about the sign of the chemical potential: 


M = 


d-2p T 
d N ' 


( 9 ) 


E = T + V. 


( 5 ) 
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The VK criterion deals with continuous families of sta¬ 
tionary localized states parametrized by their norm. The 








scaling properties of the constituents of Eq. (|^ imply that 
members of such a family will be connected by scaling 
transformations. Indeed, it is easy to show that if there 
exists a stationary solution 'ijji{r )—corresponding to the 
chemical potential fj,i and the number of atoms Ni —of 
the stationary wave equation in Eq. ([^, then 


0(r) = A(V)-P/2^i(r/A(V)) 

(10) 

= A(V)-Vi , 

(11) 


with 


= I w. 


( 12 ) 


is also a solution. According to the VK criterion, a member 


of the family of localized stationary solutions in Eqs. (10 111 
is dynamically stable—in particular against collapse or dis¬ 
persion to infinity—if its chemical potential decreases with 
the number of particles. 




(13) 


where the derivative is taken along the family. Using the 
results in Eqs. it is easy to show that in our case, 

dfi p(2p — d) Ti f ^ \ 
dN^~ d{p-d) Nf \Wi) 

So Eq. 0 supports stable stationary solitary waves if 

p > d or 0<p<^. (14) 


The second case does not currently seem physical. Recall 
that p is required to be even; then, even for d = 3, there 
is no p yielding the second inequality. 

The heuristic argument for the same conclusion goes 
as follows: let £ be the typical length scale of the atomic 
cloud. On dimensional grounds, the kinetic energy scales 
as ~ -VXjP^ while the interaction energy, which is propor¬ 
tional to the density, scales as ~ The matter wave 

will be stable if the total energy has a minimum. 

Consider first the case p > d. Then the kinetic energy 
dominates for small £(£!—>■ -foo), and the interaction 
energy for large £ {E —>• 0). Consider the energy as a 


function of £, and assume, as is reasonable, that it is a 
smooth function of it. Then dE/d£ is negative for small £, 
and positive for large £; thus, for at least one value of £, 
dE/d£ switches the sign from negative to positive, so this 
value of is a local minimum. dE/d£ could cross zero at 
multiple values of £, but at the smallest and the largest of 
these, dE/d£ must cross from negative to positive values, 
in order to match the asymptotic behavior. In any case, 
there are local minima; moreover, because E(£) is smooth, 
it follows that at least one of these local minima must be 
a global minimum. 

Now assume p < d. To show that E{£) can have no 
global or local minimum, additional assumptions will be 
needed: the reason is that while dE/d£ must in this case 
switch sign from positive to negative at least once (and 
if more than once, the first and the last switch must be 
from positive to negative, to match the asymptotic be¬ 
havior), there could, conceivably, also exist intermediate 
zeros of dE/d£ where the switch is from negative to posi¬ 
tive. These would be local minima, which would mean a 
potentially stable matter wave. However, heuristically, it 
seems unlikely that E{£) could have such an intricate struc¬ 
ture, with at least three local extrema. We can strengthen 
this argument by assuming that the functional form of 
the matter field wavefunction '0 must belong to a one- 
parametric family, parametrized by £. Then, on dimen¬ 
sional grounds, we must have 0(r) = p(vj£) j, where 
ip[z) does not depend on £. From Eqs. ([5][^ it then fol¬ 
lows that E{£) = A/£'P — B/£'^, where A and B are positive 
and do not depend on £. Then E'{£) has a single zero, at 

£q = (dB^ fQj. which 

/dB\^ 

E"i£o) = ip-d)pA^—j . 

So, A p < d, E''{£o) < 0, so £o is a global maximum, and 
there are no local or global minima. This seems the extent 
to which the heuristic argument can be pushed, and is one 
reason why it is useful to have the rigorous VK criterion. 

The standard ID nonlinear Schrodinger equation {p = 


3 






2, d = 1) is known to support solitons [TTJ [TH [H [3] and 
indeed it satisfies the first inequality in Eq. ( [T4| ). To the 
contrary, in 3D, the localized structures are known to col¬ 
lapse [Ml E]: and indeed, the combination p = 2 and 
d = 3 violates both inequalities in Eq. ( [T4| . We will now 
suggest a way to boost the dispersion sharpness to p = 4 
using dispersion management in optical lattices, and in so 
doing stabilize the 3D solitary waves. 

3. Realizing a quartic dispersion law in a shaken 

optical lattice 

Consider a single atom in a shaken 3D optical lattice: 


where the operators ~ ^ ^ = 

«:^VE(sin kx) {cos urt) are regarded as an “unperturbed Hami- 
tonian” and a “perturbation,” respectively. The corre¬ 
sponding time-dependent Schrodinger equation supports 
Floquet-type solutions, 

-1-00 \ 

, (17) 

— — OO / 

where E is the Floquet quasi-energy. The states {x)) = 
Sn^^oo the eigenstates of the Floquet Hamil¬ 

tonian: 

ni-)\^i-){x))=£^-)\x^-){x)) (18) 


ir 

H=- - \-W cos[«;(ra — ^coswt)] , (15) 

2m 

ct—x, y, z 

where m is the atomic mass, 2W the lattice depth, k = 
27r/aiat the lattice wavevector, oiat the lattice spacing, oj 
the shaking frequency, ^ the shaking amplitude (shaking 
is applied in the grand-diagonal direction), and = 

X, y, 2 . We look for the solutions of the corresponding 
time-dependent Schrodinger equation using a separation- 
of-variables ansatz: (/)(r, t) = y ^4'^°‘Kra, t). The 

time-dependent Schrodinger equations are identical for each 
of the three factors. From now on, let us concentrate on 
the time evolution of . The results can be trivially re¬ 
cast for and and we will assemble them into a 
single 3D solution in the end. 

The time-evolution equation for is 

t) = t) (16) 

at 

with 

r? 

ff{x) _ (;os[ft:(a; — ^coswt)] , 

2m 

where enters the full Hamiltonian through H = + 

jj(y)-l-7jW. Assume that ^ ^ aiat and Taylor expand 
to the first power in 

•'2 

^ + WCOSKX 
2m 

-\- (sin Kx) (cos ujt), 


with 


T-L^^) =-^ + W{cos Kx) + nhw\n){r 

2m 


+ O 0 


+ <50 


-«:^lT(sinKa;) (|n-I-l)(n|-I-|n)(n-I-1|) . (19) 


Here, the states \n) lie in a Floquet Hilbert space. 


The Floquet Hamiltonian in Eq. (191 is also periodic 
in space. Following the Bloch theorem, we will be looking 
for the eigenstates of the Floquet Hamiltonian that have 
the form 

/ +00 ^ \ 


|^(D.-ffx( 2 .)) ^ j ^ exp[i^ra] ) e' 


,iKxX 


\^— — c 


where is the cc-component of the Bloch vector. 

Denote the eigenstates and the eigenvalues of the Flo¬ 


quet Hamiltonian in Eq. (19), labeled by indexes h = 
0, ±1, ±2, ... and s = 0, 1, 2, ..., as and 

g{x),n,s ^ respectively. For zero shaking, the eigen¬ 
states and eigenenergies are related to the eigenstates 
and the eigenenergies Eq^^’^{Kx) of the sta¬ 
tionary lattice as: 


lx(“)’'^"’"'’®(a^)) ^ x‘h"'^’^"’’'{x)\n = n) (20) 


E^^^^’%K,) + hu;n . 


( 21 ) 


For a shaking frequency close to the transition frequency 
between the ground and the first excited energy bands in 


4 




the middle of the Brillouin zone. 


1 




( 22 ) 


these bands hybridize [isiiiaiiTiiiiiii]. Observe that for 
a deep blue detuning, the bands will exchange roles and the 
parabolic dispersion law in the ground band will become 
an inverted parabola. By continuity, for each shaking am¬ 
plitude, there will exist a shaking frequency for which the 
quadratic term in the dispersion law vanishes, the quar- 
tic term becoming dominant. Below, we analyze a sample 
set of parameters where this phenomenon indeed happens. 


The resonance condition in Eq. (22) allows us to use de¬ 
generate perturbation theory, truncating the Hilbert space 
to a 2D space spanned by the states = 

and = —1), with the unperturbed energies 

j^(x), (rf)^ j j^(x), s=i )—huj, respectively. Our con¬ 

struction of a 3D quartic dispersion curve is completed by 
the observation that the identical dispersion law will be 
generated in the other two spatial directions, and that the 
resulting dispersion hypersurface is simply the sum of the 
three. 


4. An outline of a scheme for experimental real¬ 
ization 

Let us now discuss a sample experimental realization 
of the above scheme. Consider a set of parameters in¬ 
spired by the scheme used to induce a long-range ferro¬ 
magnetic order in an atomic gas [19], as follows. N = 10^ 
^^^Cs atoms are loaded in a 3D optical lattice with lat¬ 
tice spacing oiat = 532 nm and depth 2W = 6.8 Er, where 
Er = /2m is the recoil energy, and k = n/2 = 7r/aiat 

is the wavevector of the lattice light. The interactions are 
repulsive and some trapping will also be required at this 
stage. We are going to postpone the choice of the scatter¬ 
ing length till later. Initially, there is no shaking. Then 
a low-frequency shaking starts, with the amplitude 
in the grand diagonal direction, where f = 0.0056 oiat is 
the shaking amplitude along each of the three ID-lattice 


direction. Subsequently, the shaking frequency is slowly 
ramped to w = 4.75 Er//!, a frequency slightly to the red 
from the interband spacing in the center of the Brillouin 
zone. For the same frequency, the edges of the zone un¬ 
dergo a complete adiabatic population inversion during 
the ramp. The final shaking frequency is tuned in such 
a way that the quadratic terms in the dispersion law of 
the—heavily hybridized—ground band are canceled: the 
remaining curve is the 3D inverted quartic parabola de¬ 
picted in Fig. Following the analogy with the ID gap 
solitons [20], we concentrate on the time-evolution of the 
complex conjugate of the wavefunction. Its time evolution 
is governed by a non-inverted quartic dispersion law: 


at 

77 > 0, g < 0, 


q4^ q4^ 


+ 


+ 


dx^ dy'^ dz^ 


(23) 


where tk = 4/*, Tj = —rj^ and g = —g. Accordingly, the 
stability of the solitary wave will require repulsive inter¬ 
actions; as an example, we set the scattering length to 
as = 2.4aBohr- At this point, trapping can be removed, 
following a numerically optimized schedule, and the re¬ 
sulting atomic cloud should be entirely self-supporting. 
The dispersion constant entering the time-dependent wave 


equation in Eq. (231 will have the value 


7 ) = -4! 


while the coupling constant in Eq. (23) will be given by 
g = —X (ATTh? as/m), where 

/ + aiat/2 
-aiat/2 

here a is any of the three Cartesian directions, cc, y, or 
z. Using our chosen set of parameters, this works out to 


fj = 0.3ER(aiat/S) and g = -0.6 (see Eqs. (201 


and (21) for notation). For the above set of parameters, a 


Gaussian variational estimate predicts the wavevector dis¬ 
tribution diameter of AK = 2((A^))^/"^ « 0.018(27r/aiat); 
for this degree of momentum localization, the solitary wave 
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Figure 1: (color online). Floquet energies (red dashed 

lines) and ~^’^{Kx) (blue dotted line), the adiabatic descen¬ 

dants of the ground and the first excited bands, respectively, as func¬ 
tions of the Bloch vector Kx- The parameters used correspond to a 
system of ^^^Cs atoms in an optical lattice with depth 2W = Q.SEji 
and spacing aiat = 532 nm, shaken in the grand-diagonal direction 
with frequency cj = A.75 Ej^/h and amplitude y/S x 0.0056 d. The 
curves are calculated using degenerate perturbation theory applied 
to the Floquet Hamiltonian. For our set of parameters, quadratic 
terms in the ground band dispersion are canceled, and the formerly 
sub-leading inverted quartic dispersion (the red solid line in the inset) 
becomes dominant. Dispersion laws in the two remaining directions, 
y and 2 , are identical to the one along x. Inset: a magnified view of 
the peak of the Floquet energy ^{Kx) (red dashed line), show¬ 

ing that it is very close to the pure—formerly sub-leading—inverted 
quartic dispersion (the red solid line). The axes labels of the main 
plot are also the axes labels for the plot in the inset. 

would occupy only ~ l/20th of the Brillouin zone in each 
of the three directions, and thus it is expected to be fully 
mobile. (Above, we used the exact result in Eq. ^ to re¬ 
late the mean fourth power of the quasi-momentum to the 
total energy.) This set of parameters also makes the inter¬ 
atomic dipole-dipole interaction small: its relative strength 
compared to the contact interaction is Cdd = -25 (see p. 70 
of Ref. HD). 


5. Conclusion and outlook 

We have shown that a 3D Bose-Einstein condensate 
with a quartic dispersion law satisfies the VK criterion for 
the stability of localized structures in a dispersive non- 
relativistic self-focussing medium, and that such a disper¬ 
sion law can be realized in shaken optical lattices umia 
imiiiiis]. Experimentally, as in the realization of band- 
gap solitons EH], our scheme will involve an inverse quartic 
dispersion and repulsive condensates. 

We see rotation sensing as the primary area of ap¬ 
plication of 3D solitary waves. Interferometry with self- 
attracting waves was a success from the onset [22] (see 
also [231 EH EH Uni)) showing a manyfold increase in the 
fringe visibility as compared to free waves. Further im¬ 
provements in sensitivity are expected for slow splitting 
of solitary waves on barriers |27j where the appearance of 
macroscopic quantum superpositions is predicted [531ES] ■ 
But while linear solitary wave interferometers are ideally 
suited for acceleration sensing, Sagnac measurements re¬ 
quire a multi-dimensional geometry. There are waveguide- 
based proposals [30]; their principal drawback—in addi¬ 
tion to the fact that the loading of atoms into the guide 
is difficult, and that the requirements for the stability of 
the guiding fields are very strict—is the inflexibility of 
the interferometric scheme. In developing an interfero¬ 
metric scheme, one starts with a seed idea and then builds 
upon it by successive modifications-usually strengthening 
the beamsplitters in some way—that result in increases 
in the interferometric area (space-time for the gravitome- 
ters/accelerometers, and space-space for the rotation sen¬ 
sors). For waveguide-based schemes, this requires rebuild¬ 
ing the waveguide each time, with all the accompanying 
optimization and stabilization issues. In contrast, schemes 
based on free-space propagation allow for easy implemen¬ 
tation of improvements. The present work paves the way 
for such schemes, by proposing a method to obtain solitary 
waves that propagate in free space. 
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